On the excitation of f-modes and torsional modes by magnetar giant 
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ABSTRACT 

Magnetar giant flares may excite vibrational modes of neutron stars. Here we compute an estimate of initial 
post-flare amplitudes of both the torsional modes in the magnetar's crust and of the global f-modes. We show 
that while the torsional crustal modes can be strongly excited, only a small fraction of the flare's energy is 
converted directly into the lowest-order f-modes. For a conventional model of a magnetar, with the external 
magnetic field of ~ 10 15 G, the gravitational-wave detection of these f-modes with advanced LIGO is unlikely. 



Subject headings: Neutron stars 



Introduction 



The gamma- and x-ray flares from Soft Gamma Re- 
peaters (SGRs; Mazetz et al. 1979, Hurley et al. 1998, 
2004) are believed to be powered by a sudden release of 
magnetic energy stored in their host magnetars (Thomp- 
son & Duncan 1995). An SGR flare may excite vibrational 
modes of a magnetar (Duncan 1998). Indeed, torsional 
oscillations of a magnetar provide an attractive explana- 
tion for some of the quasi-periodic oscillations (QPOs) ob- 
served in the tails of giant flares (Barat et al. 1983, Is- 
rael et al. 2005, Strohmayer & Watts 2005, van Hoven 
& Levin 2011, Gabler et al. 2011, Colaiuda & Kokkotas 
2011). Moreover, there is now some evidence for QPOs 
during the normal, non-giant flares in SGR 1806-20 (El- 
Mezeini & Ibrahim 2010). 

Excitation of low-order f-modes is also of considerable 
interest, because of the f-modes' strong coupling to po- 
tentially detectable gravitational radiation. The sensitiv- 
ity of the ground-based gravitational-wave interferometers 
has dramatically improved over the last 5 years (Abott 
et al. 2009a, Acernese et al. 2008), and interesting up- 
per limits on the f-mode gravitational- wave emission from 
the 2004 SGR 1806-20 giant flare, a possible 2009 SGR 
1550-5418 giant flare, and several less energetic bursts 
have recently been obtained (Abott et al. 2008, Abott et 
al. 2009b, Abadie et al. 2010, see also Kalmus et al. 2009). 



Advanced LIGO and VIRGO are expected to become op- 
erational in the next 5-7 years, and it is of interest to 
predict the strength of expected gravitational-wave signal 
from future giant flares. In this paper, we compute a the- 
oretical estimate for the amplitude of the torsional and 
f-modes expected to be excited in a giant flare. We show 
that only a small fraction of the flare energy is expected 
to be pumped into the low-order f-modes, and estimate 
the signal-to-noise ratio for the future giant flare detection 
with advanced LIGO. By contrast, the torsional modes can 
be strongly excited and may well be responsible for some 
of the observed QPO's in magnetar flares. 

2. The general formalism 

The giant flares release a significant fraction of the free 
magnetic energy stored in their host magnetars. Two dis- 
tinct mechanisms for this have been proposed: 1. Large- 
scale rearrangement of the internal field, facilitated by a 
major rupture of the crust (Thompson & Duncan 1995, 
2001; we shall refer to it as the internal mechanism, IM), 
and 2. Large-scale rearrangement of the magnetospheric 
field, facilitated by fast reconnection (Lyutikov 2006, Gill 
& Heyl 2010; we shall refer to it as the external mechanism, 
EM ). Both processes may well be at play: the IM would 
likely serve as a trigger for the EM (however, as was argued 
in Lyutikov 2003, EM may also be triggered by slow motion 
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of the footpoints of a magnetospheric flux tube, leading to 
a sudden loss of magnetostatic equilibrium). Observation- 
ally, the extremely short, a few microseconds rise time of 
the 2004 giant flare in SGR 1806-20 (Hurley et al. 2004) 
gives reason to believe that EM was at play in that source: 
the IM operates on a much longer Alfven crossing timescale 
of 0.05-0.1 seconds. The long timescale for the IM implies 
that it would not be efficient in exciting the f-modes which 
have frequencies of over a kHz; this was recently indepen- 
dently emphasized by Kashiyama & Ioka (2011). We shall 
come back to this point more quantitatively at the end 
of this section. In what follows we will thus concentrate 
on the mode excitation by the EM. The IM may well be 
efficient in exciting lower-frequency torsional modes; how- 
ever, as we discuss below, its coupling to these modes is of 
similar strength to that of the EM, and as we will show, 
both mechanisms may be efficient in exciting the large- 
amplitude torsional modes. 

During the large-scale EM event, the magnetic stresses 
at the stellar surface change instantaneously (in several 
light-crossing timescales) by, at mosl0, order 1. We shall 
characterize the change of the magnetic stress by the 3 
components 
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where B is some characteristic value of the surface mag- 
netic field and f r , and are functions are of order 
1 in the strongest possible flares and are smaller for the 
weaker flares. Consider now a normal mode of the star 
with an eigenfrequency Lo n and a displacement wavefunc- 
tion t; n (r,6,4>). We treat the changing surface magnetic 
stress as an external perturbation acting on the mode. The 
Lagrangian of the free (pre-perturbation) mode is given by 



Lf rC c( a m On) — ~^ m nQ> n 



(2) 



where a n is the generalized coordinate corresponding to 
the normal mode, m n is the effective mass given by 



d A rp{r)&{r) 



and p(r) is the density. The Lagrangian term characteris- 
ing the mode's interaction with external stress is given by 



1 During the short timescale of the EM event no new flux can emerge 
from the star. 



(cf. section 2 of Levin 1998) 



where 
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and the displacement £ is evaluated at the radius of the 
star R. The full Lagrangian for the nth mode is given b>H 



_L(d n , (i n ) — Z/f rcc -\- F ma g(x n . 



where 
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is the characteristic energy stored in the star's magnetic 
field and and a n is the coupling coefficient given by 



n = / C (R, 0, 0) • f(9, 4>) sin Offldt, 



where 



(8) 



/ = f r (9, </>)Sr + f e (0, 0)e e + U(9, fie,/, (9) 



It is now trivial to find the motion resulting from the sud- 
den introduction of the external stress at moment t = 0. 
The coordinate a„ oscillates as follows: 



a n (t) = a n [1 - cos (<jJ n t)\ , 
where the amplitude is given by 

_ ^n-S'mag 
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The energy in the excited mode is given by 

„2 T?2 
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We now briefly revisit the mode excitation by the IM. 
In this case, the interaction Lagrangian of a mode with 
the magnetic field is described by the following volume 
integral: 

L int =a n ( d 3 rf L (r)-Z n {r), (13) 



(3) where f L 



V x B 



x B is the lorentz force per unit vol- 



ume. Since Jl ~ B 2 /R, one can see that the coupling of 



2 We work in the linear regime, and don't take into accounf the non- 
linear coupling between the modes. The mode amplitudes <C 1 found 
at the end of our calculation indicate that this is a good approxima- 
tion. 
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the internal field variation to the mode is of the same order 
of magnitude as that of the external field variation, pro- 
vided that the external and internal fields are of the same 
order of magnitude. However, the IM mechanism acts on 
a much longer timescal^l TAif VO n ~ 0.1s than the typical 
f-mode period of tj ~ 0.0005s, so the f-mode oscillator 
would be adiabatically displaced without excitation of the 
periodic oscillations. One can show that the typical sup- 
pression factor of the IM relative to the EM excitation is at 
least of order 27TTAif ve n/V/ in the mode amplitud^]. This 
factor is so large that even if internal field was stronger 
than the external field by an order of magnitude, the IM 
excitation would still be suppresed relative to the EM one. 

Is there a way around this suppression factor? Poten- 
tially, IM could feature a collection of many localized MHD 
excitation, with the timescale for each one being deter- 
mined by the Alfven-crossing time of each of the exci- 
tation domain. If the domains were small enough, their 
timescales could be more closely matched with the f-mode 
period (Melatos, private communications). However, in 
this case the magnitude of the overlap integral in Eq. (fT3"|) 
would be reduced by a factor ~ (R/AR) 3 , where AR is 
the characteristic size of the excited domain. The domains 
would contribute incoherently to the amplitude of the ex- 
cited mode, thus the contribution of an individual domain 
would have to be multiplied by (R/ ' AR) 3 / 2 in the (some- 
what unlikely) limit where the active domains occupy the 
whole star. Thus, while the timescale of the mini-flares 
could be well-matched with the f-mode period, their over- 
all conrtribution to the overlap integral in Eq. (|13p would 
be suppressed by ~ (R/ AR) 3 / 2 . In the optimal case that 
the mini-flares have the same timescale as the f-mode pe- 
riod, R/ AR ~ TAifvcn/17. Therefore, the collection of 
mini-flares would not give us any gain in the mode ex- 
citation amplitude, as compared to the IM estimate given 
in the previous paragraph. 

Two applications of the formalism for the mode excita- 
tion by the EM mechanism developed above are presented 



3 This timescale could be shorter if the supcrfiuid neutrons are de- 
coupled from the MHD (Easson & Pcthick 1979, van Hoven & Levin 
2008, Andersson, Glampedakis. & Samuelsson 2009). However, given 
the large fluid velocities expected in IM, this decoupling seems un- 
likely. 

4 This can be formalized by the following argument: consider a har- 
monic oscillator of proper frequency u>o, initially at rest, which is 
externally driven by force f(t). The amplitude of the induced oscil- 
lation at the proper frequency is proportional to /(i^o), the Fourier 
transform of f(t) evaluated at u>q. For a step function, as expected in 
EM, f(uj) oc l/ui. On the other hand, for a smooth pulse of duration 
t, as expected in IM, the Fourier transform is suppressed and scales 
at most as f(u>) tx (wt) _1 1/w when ujt S> 1. 



in the next two sections. 

3. f-modes and gravitational waves 

In order to estimate an effective f-mode mass, we have 
computed the I = 2 f-mode displacement functions for a 
neutron staj^l in the Cowling approximation. Convenient 
scalings are 

m n = q M M, 
2 GM 

£rCR,M) = a n Y 2m (e,<t>). 

(14) 

In our fidutial model qm = 0.046, where we have nor- 
malised the mode wavefunction so that e r ■ £,2m(R, 6, 4>) = 
Y2 m (9,4>). Our reference number — 1.35 was obtained 
using a fitting formula for fully relativistic f-mode frequen- 
cies^! from Andersson & Kokkotas (1996). The amplitude 
of the f-mode is given by 

where 

GM 2 

-kgrav = — (16) 

is of the same order as the gravitational binding en- 
ergy of the neutron star. For our fidutial model with 
B = 10 15 G, M = 1.4M , and R = 10km, we get 
0,2m /R ~3x 10~ 6 a2m- The energy in the f-mode is 

Ef = ^~^h ~ L5 x ^- 6 a 2 m E mas (17) 

^HmQuj -^grav 

for our fidutial parameters. This energy is drained from 
the star primarily through emission of gravitational waves. 



5 We constructed our neutron star model using the equation of state 
from Douchin & Haensel (2001) and Haensel & Pichon (1994). In 
calculating the f-mode we treated the whole star as a fluid, neglecting 
the effects of bulk- and shear moduli. 

6 We are not being consistent in, on the one hand, using the Cowling 
approximation for a Newtonian star to determine the effective mode 
mass, but on the other hand using the published relativistic calcu- 
lations for the mode frequencies. Normally, Newtonian calculations 
would be sufficient, given the many unknown details of the flare and 
the many poorly constrained parameters we'd already introduced into 
the model, and the formalism we developed in the previous section is 
manifestly Newtonian (but can be generalized to relativistic regime if 
the need arises). However, as we show below, the signal-to- noise ra- 
tio for the gravitational-wave detection is very sensitive to the mode 
frequency, and therefore we try to be accurate in characterizing these 
frequencies. 
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The total amount of energy carried by gravitational waves 4. Torsional modes 



is therefore 
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where / = w„/27r is the f-mode frequency in Hz, {h 2 } 
is the direction and polarisation averaged value of the 
square of the gravitational-wave strain h as measured by 
observers at distance d from the source. This expres- 
sion allows us to estimate the expected signal-to-noise ra- 
tio for ground based gravitational wave interferometers 
(cf. Abadie et al 2010). One can use the fact that nearly 
all the gravitational-wave signal is expected to arrive in 
a narrow-band around the f-mode frequency, and that the 
signal form (the exponentially-decaying sinusoid) is known. 
The Wiener-filter expression for the signal-to-noise can be 
written as 
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(20) 



where h(f) is the Fourier transform of the time-dependent 
gravitational wave strain h(t). As is standard for narrow- 
band signal, we have used Parseval's theorem to convert 
the integral over / to the integral over t from Eq. (|18l) . 
and, following Abadie et al. 2010, we have approiximated 
h 2 with the average (h 2 ). At frequencies of a few kHz the 
spectral density, Sh(f), of the ground based detectors like 
Advanced LIGO and Virgo is dominated by shot-noise and 
is proportional to f 2 . This makes the signal-to- noise ratio 
for observations of magnetar f-modes excited in giant flares 
particularly sensitive to frequency (oc /~ 3 ). For Advanced 
LIGO we find 



S nn ~ /2000 Hz 
N * °-° 7 a2m ( — — 



D 



10 15 G 

2 



I? 



10 km 



1 kpc 
d 

0.07 Mq 



x 

1/2 



(21) 



Here we used tabulated^ Sh(f) from the LIGO document 
ligo-t0900288, which gives 

S h (f) = 8.4- lO-^Hz" 1 (//2000 Hz) 2 for the shot-noise 
dominated part of the curve. 



7 These sensitivity curves represent the incoherent sum of principal 
sources of noise as they are currently understood. 



Intuitively, one expects torsional modes to be strongly 
excited during the magnetar flares (Duncan 1998), since it 
is the free energy of the twisted magnetic field that is being 
released. These have much lower proper frequencies than 
the f-modes (with the fundamental believed to be in the 

range 10 40Hz, see Steiner & Watts 2009 and references 

therein), which can be well-matched to the Alfven frequen- 
cies inside the star. Thus both EM and IM are likely to 
play a role in the torsional mode excitation. Here, we con- 
sider the EM explicitly but keep in mind that IM would 
give a similar answer. 

For the torsional modes in the crust, the displacement 
is given by 



(, n im{r, 0, 4>) = g n {r)r x VYi 



(22) 



and it is convenient to normalise the wavefunctions so that 
9n(R) — 1- Here n — 0, 1, ... is the number of radial nodes. 

With this normaliation, the effective mode mass m n i m ~ 
TOcrust ~ 0.01M, and from Eq. (fTT|) one gets for the mode 
amplitude normalized by the star radius: 



a n i 



( B \ 2 R 0.014M Q /100Hz\ 2 



m ni i. 



(23) 

Thus we see that for a reasonable range of parameters it is 
feasible that the crustal torsional modes would be strongly 
excited by a giant flare. 

4.1. Magnetic modes 

Recently, Kashiyama & Ioka (KI, 2011) suggested that 
certain types of MHD modes that may be strongly excited 
during a giant flare, are coupled to gravitational radia- 
tion and may therefore become an interesting source for 
advanced LIGO. KI focus on the polar modes of Sotani 
& Kokkotas (2009); the MHD modes found by Lander & 
Jones (2011a,b) also satisfy some of the KI's criteria. 

While interesting, this idea has potential caveats that 
need further investigation. The polar, axially-symmetric 
modes of Sotani & Kokkotas (2009) exist in idealized poly- 
tropic stars; in a real neutron star the frequency of such 
modes would likely be many times greater due to the 
strong stable stratification (Goldreich & Reisenegger 1992, 
Reisenegger & Goldreich 1992) and their nature would 
likely be closer to that of the g-modes. More importantly, 
KI assume that the oscillations are long-lived, ~ 10 7 os- 
cillation periods. However, MHD modes are notoriously 
capricious. While the idealized modes of Sotani & Kokko- 
tas (2009) and Lander & Jones (201 la, b) are protected by 
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symmetry, the global magnetic modes in more realistic con- 
figurations may couple to a variety of localized Alfven-type 
modes, and may thus be quickly damped via phase mix- 
ing and resonant absorption (Goedbloed & Poedts 2004, 
van Hoven & Levin 20f 1 ). Thus, in our view, there is cur- 
rently no compelling reason to believe that the magnetic 
modes can be substantially longer lived than the observed 
magnetar QPOs. 

5. Discussion 

In this paper, we have computed the excitation of the 
f-modes and crustal tosional neutron-star modes by a giant 
flare. Corsi & Owen (2011) recently computed the mag- 
netic energy that can be released during the flar^f], and 
found values comparable to E mag . However, in this work 
we showed that only a small fraction of the released flare 
energy is converted into the f-modes, and that the associ- 
ated gravitational- wave emission is correspondingly weaker 
than has been previously hoped (cf. Abadie et al. 2010 and 
Corsi & Owen 2011). From Eq. (|2"Tj) . our fiducial model 
does not look promising for future advanced LIGO detec- 
tion of a giant flare, even if a.2 m ~ 1, i.e. if the released 
electro-magnetic energy is of order of the total magnetic 
energy of the star, ~ 10 47 erg (the most energetic of the 
3 observed giant flares released few x 10 erg). However, if 
the surface field is significantly larger than 10 15 G and/or 
the star is greater than 10km in radius (which would re- 
duce the f-mode frequency and increase the contact surface 
area), then one can become more hopeful about the poten- 
tial detection. 

On the other hand, we have seen that there is no diffi- 
culty in exciting the crustal torsional modes to a large am- 
plitude. Whether or not this leads to the observed quasi- 
periodic oscillations in the flare's tail (Israel et al. 2005, 
Strohmayer & Watts 2005, Watts & Strohmayer 2006) de- 
pends crucially on the dynamics of hydromagnetic cou- 
pling between the crustal modes and the Alfven modes of 
the magnetar core (Levin 2006, 2007, van Hoven & Levin 
2011, Gabler et al. 2011, Colaiuda & Kokkotas 2011). 
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